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' In this article we present deep, high-resolution radio interferometric observations at 153 MHz to complement the 

I extensively studied NOAO Bootes field. We provide a description of the observations, data reduction and source catalog 

• , construction. From our single pointing GMRT observation of ~ 12 hours we obtain a high-resolution (26" x 22") image 

' of ~ 11.3 square degrees, fully covering the Bootes field region and beyond. The image has a central noise level of 

~ 1.0 mjy beam ^, which rises to 2.0-2.5 mjy beam ^ at the field edge, placing it amongst the deepest ~ 150 MHz 
surveys to date. The catalog of 598 extracted sources is estimated to be ~ 92 percent complete for > 10 mJy sources, 
while the estimated contamination with false detections is < 1 percent. The low RMS position uncertainty of 1.24" 
facilitates accurate matching against catalogs at optical, infrared and other wavelengths. Differential source counts are 
^ ' determined down to < 10 mJy. There is no evidence for flattening of the counts towards lower flux densities as observed 

in deep radio surveys at higher frequencies, suggesting that our catalog is dominated by the classical radio-loud AGN 
population that explains the counts at higher flux densities. Combination with available deep 1.4 GHz observations 
yields an accurate determination of spectral indices for 417 sources down to the lowest 153 MHz fiux densities, of which 
^ ' 16 have ultra-steep spectra with spectral indices below —1.3. We confirm that flattening of the median spectral index 

I towards low fiux densities also occurs at this frequency. The detection fraction of the radio sources in NIR Ks-hand is 

, found to drop with radio spectral index, which is in agreement with the known correlation between spectral index and 

On ' redshift for brighter radio sources. 
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1. Introduction pear to hold down to very low flux levels (e.g., Afonso et al. 

2011). Concentrating on the faintest sources from surveys 
Surveying the radio sky at low frequencies (< 300 MHz) is a ^^adc at the lowest frequencies is therefore an obvious way 
unique tool for mvestigatmg many questions related to the pushing the distance limit for HzRGs beyond the present 
formation and evolution of massive galaxies, quasars and highest redshift of TN J0924-2201 at z = 5.1 (van Breugel 
^ . clusters of galaxies (e.g., Miley & De Breuck 2008). Low- al. 1999) and probing massive galaxy formation into the 
■ frequency radio observations benefit from the steepness of epoch of reionization 
radio spectra of various types of cosmic radio sources, such Galaxy clusters containing diffuse radio sources appear 
as massive HzRGs (high-redshift radio galaxies; redshift ^ ^g^^^ i^^^^ X-ray luminosities and galaxy velocity disper- 
z>2) and^diffuse halo & relic emission m nearby galaxy ^-^^^ (g_g^^ Hanisch 1982), which are thought to be char- 
clusters [z ^ 0.1). acteristics of cluster merger activity (e.g., Giovannini & 
HzRGs arc amongst the most massive galaxies in the Feretti 2000; Kempner & Sarazin 2001). Synchrotron halos 
early Universe (e.g., Miley & De Breuck 2008), usually lo- and relics provide unique diagnostics for studying the mag- 
cated in forming galaxy clusters with total masses of more netic field, plasma distribution and gas motions within clus- 
than 10^"^ Mq (e.g., Venemans et al. 2007). The most ef- ters, important inputs to models of cluster evolution (e.g., 
ficient way of finding HzRGs is to focus on USS (ultra- Feretti & Johnston-HoUitt 2004). Cluster synchrotron emis- 
steep spectrum; S,, oc with a < -1) radio sources sion is known to be related to the X-ray gas and pinpoints 
(Roettgering et al. 1997; De Breuck et al. 2002). This was shocks in the gas. Further, cluster radio emission usually 
reinforced by Klamer et al. (2006) who showed that the has steep radio spectra {a < -1), the radiating electrons 
radio spectra of HzRGs in general do not show spectral are old and can provide fossil records of the cluster history 
curvature, but are straight. The USS selection criteria ap- (e.g., Miley 1980). 

Several low-frequency surveys have been performed in 
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et al. 1959; Bennett 1962; Pilkington & Scott 1965; Gowcr 
et al. 1967; Hales ct al. 1988, 2007), the UTR-2 sky sur- 
vey between 10-25 MHz (Braude ct al. 2002), the VLSS 
at 74 MHz (Cohen et al. 2007) and the ongoing MRT 
sky survey at 151.5 MHz (e.g., Pandey & Shankar 2007). 
These surveys are limited in sensitivity and angular res- 
olution, mainly due to man-made RFI, ionospheric phase 
distortions and wide-field imaging problems. Recent devel- 
opments in data reduction techniques make it possible to 
perform deeper surveys (< 50 mJy bcam~^) of the low- 
frequency sky at higher resolution (< 30"). One telescope 
that significantly improved this situation is the GMRT. A 
few deep, single-pointing surveys at 153 MHz (e.g., Ishwara- 
Chandra & Marathe 2007; Sirothia et al. 2009; Ishwara- 
Chandra et al. 2010) have been performed, yielding noise 
levels between 0.7-2 mJy beam^^ and resolutions between 
20-30". And the TCSS^ is a new, ongoing 153 MHz GMRT 
sky survey, aimed at covering the full northern sky down to 
DEC > -30° at a 20" resolution and a 7-9 mJy beam^^ 
noise. 

In this article we present deep, high-resolution GMRT 
observations at 153 MHz of the NOAO Bootes field. The 
Bootes field is a large (~ 9 square degree) northern field 
that has been targeted by surveys spanning the entire 
electromagnetic spectrum. This field has been extensively 
surveyed with radio telescopes including the WSRT at 
1.4 GHz (de Vries et al. 2002), the VLA at 325 MHz 
(Croft et al. 2008), and will be complemented with deep 
74 MHz EVLA observations. The large northern NDWFS 
survey (Jannuzi & Dey 1999) provided 6 colour images 
[By^RIJHK) to very faint optical and NIR flux lim- 
its. Additional, deeper NIR images in J Kg are available 
from the FLAMEX survey (Elston et al. 2006), while addi- 
tional z-band images are available from the zBootes cam- 
paign (Cool 2007). The entire area has also been surveyed 
by Spitzer in seven IR bands ranging from 3.6 to 160 /im 
(Eisenhardt et al. 2004; Houck et al. 2005). Chandra has 
covered this area in the energy range of 0.5-7 keV to a 
depth of ~ 10"^'' ergs cm~^, yielding 3200 quasars and 
30 luminous X-ray clusters up to rcdshift z ^ 1 (Murray 
et al. 2005; Kcnter et al. 2005). The UV space telescope 
GALEX has covered the Bootes field. AU the 10,000 galax- 
ies brighter than R < 19.2 and X-ray / IRAC / MIPS 
QSOs brighter than R < 21.5 have redshifts through the 
AGES project (Kochanek et al. in preparation). Based on 
the shallow Spitzer data, 3 HzRGs with photometric red- 
shifts of z > 4 have been identified in the Bootes field (Croft 
et al. 2008). Also, Cool et al. (2006) report the discovery of 
3 quasars with spectroscopic redshifts z > 5, while McGreer 
et al. (2006) found a quasar at z = 6.1. 

Given the size of the GMRT FoV (field-of-view - 3.5 de- 
grees) and angular resolution ('^ 25") at 153 MHz, the 
Bootes field is a well-matched region for conducting a deep 
survey. Combined with the existing multi-wavelength sur- 
veys, our deep 153 MHz Bootes field observations allow for 
a complete study of faint (> 20 mJy) low-frequency radio 
sources. For the data reduction, we used the recently de- 
veloped SPAM calibration software that solves for spatially 
variant ionospheric phase rotations (Intema et al. 2009) . In 
our initial analysis of the 153 MHz source catalog we focus 
on determining source counts down to the detection limit, 
and identifying steep-spectrum radio sources that are can- 
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Date June 3, 2005 June 4, 2005 



LST range 


12-20 hours 10-19 hours 


Local time range" 


20:00-04:00 18:00-03:00 


Time on target 


359 min 397 min 


Primary calibrator 


3C 286 3C 286 


Time on calibrator 


94 min 108 min 


Integration time 


16.8 sec 


Polarizations 


RR, LL 


No. of channels 


128 


Channel width 


62.5 KHz 


Total bandwidth 


8.0 MHz 


Target RA,DEC 


14^^32^05. 75^ +34°16'47.5" 



Indian Standard Time (1ST) — Universal Time 



(UT) + 5:30 

Table 1. Overview of GMRT observations on the Bootes 
field. 



didate HzRGs. A more extended analysis of our source cat- 
alog in combination with the other multi-wavelength cata- 
logs is planned in a subsequent article. 

In Section 2, we describe the GMRT 153 MHz obser- 
vations and data reduction. In Section 3, we present de- 
tails on the source extraction and catalog construction. 
Section 4 contains the initial analysis and discussion of the 
source population. Conclusions are presented in Section 5. 
Throughout this article, source positions are given in epoch 
J2000 coordinates. 

2. Observations and Data Reduction 

In this section, we describe the observations and data re- 
duction steps that led to the production of the 153 MHz 
image that is the basis of the survey. 

2.1. Observations 

We used the GMRT (e.g., Nityananda 2003) at 153 MHz to 
observe the Bootes field, with the observing details given in 
Table 1. There were typically 27 out of 30 antennas avail- 
able during each observing run. We used 3C 286 as fiux, 
bandpass and phase calibrator. Observation cycles of ~ 50 
minutes were split into ^ 38 minutes on the target field 
and ~ 10 scans on the calibrator. The relatively high over- 
head in calibrator observations was justified by the need 
to monitor the GMRT system stability, RFI conditions and 
ionospheric conditions, and to ensure consistency of the flux 
scale over time. 

2.2. Data Reduction 

The data reduction was performed in two stages. The flrst 
stage consisted of 'traditional' calibration, in which the flux 
scale, bandpass shapes and phase offsets were determined 
from the calibrator observations, transferred to the target 
field data, after which the target field was self-calibrated 
and imaged for several iterations. During the second stage 
of the data reduction, we made use of the recent SPAM 
software package (Intema et al. 2009) that incorporates 
direction-dependent ionospheric phase calibration. 

For the traditional calibration we used the Astronomical 
Image Processing Sofware (AIPS; e.g.. Bridle & Greisen 
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Field diameter'* 


6.8° 


Pixel size 


4" 


Weighting^ 


robust 0.5 


VV liLiU. iliidtiiiigj 


7^nl'\;' TiPn TTiTi f in fof - Qpri 1 




mn If 1- T7"pmiP'np'\7' '^■\7'Tif np'^i'^ 

lllLll ul J.J. v-vJ Llv-lH^ y L> y lHjliColD 


IMiiTnViPr of fpirpf^ 


199 


Facet diameter 


0.57° 


Facet separation 


0.47° 


Deconvolution*^ 


Cotton-Schwab CLEAN 


CLEAN box threshold 


5a 


CLEAN depth 


2a 


Restoring beam 


26" X 22" (PA 78°) 


SPAM-' calibration cycles 


2 


Peeled sources* 


24 


Layer heights (and weights) 


100 km (0.25) 




200 km (0.50) 




400 km (0.25) 


Parameter 7 (^) 


5/3 


Model parameters 


20 


Model fit phase RMS^ 


19.2° ±3.0° 


Peeling corrections applied 


yes 


Reference catalog' 


NVSS 


" We map more than twice the HPBW diameter to allow for 



deconvolution of nearby bright sources. 

Briggs (1995) 
"Perley (1989); Cornwell & Perley (1992) 

Conway et al. (1990) 
■= Schwab (1984); Cotton (1999); Cornwell et al. (1999) 
f Intema et al. (2009) 

^ Specified for the final (second) calibration cycle only. 

Power-law slope of the assumed phase structure function. 
' Condon et al. (1994, 1998) 

Table 2. Overview of wide-field imaging parameters (upper 
panel) and SPAM parameters (lower panel). 

1994) package. Initial flagging of bad baselines and ex- 
cessive visibility amplitudes (mostly RFI) resulted in a 
6.75 MHz effective bandwidth. This was followed by am- 
plitude, phase and bandpass calibration on 3C 286. We 
adopted a Stokes I flux density of 31.01 Jy on the Per ley- 
Taylor 1999.2 scalc^, which is derived from the Baars et al. 
(1977) scale for 3C 295 (the flux scale is discussed in 
more detail in Section 3.3). To reduce the data volume, 
each 4 consecutive frequency channels were combined to 
form 27 channels of 0.25 MHz each. After initial imag- 
ing of the target field (sec Table 2), the calibration on 
the target field data was improved by three rounds of 
phase-only self-calibration and one round of amplitude & 
phase self-calibration. To remove excessive visibilities at 
a lower level, the obtained source model was temporar- 
ily subtracted from the visibilities, after which the residual 
visibilities were manually and semi-automatically flagged 
per baseline. After re-imaging, the noise (RMS of the im- 
age background) in the inner half of the (uncorrected) 
primary beam area, was approximately 1.4 mJy beam~^. 
Near the brightest three sources (with apparent flux densi- 
ties larger than 1 Jy), the local noise was measured to be 
> 2.2 mJy beam~^. 

Despite the good overall quality of the self-calibrated 
image from the traditional calibration, there were sig- 

^ Defined in the VLA Calibrator Manual; http://www.vla.nrao. 
edu/astro/ca lib/manual/ baars. htm I 



nificant artifacts present in the background near bright 
sources, limiting the local dynamic range to a few hun- 
dred. To suppress artifacts related to direction-dependent 
calibration errors, we applied the SPAM algorithm on the 
self-calibrated(!) visibility data. The used SPAM parame- 
ters are reported in Table 2 (for details see Intema et al. 
2009). This unconventional combination of self-calibration 
and SPAM is used to bypass the problem of antenna-based 
phase discontinuities, as observed in the calibration so- 
lutions for several antennas. Applying SPAM after self- 
calibration can be considered as a perturbation to the over- 
all gradient correction found by self-calibration. The re- 
sulting noise in the inner half of the (uncorrected) pri- 
mary beam area is on average 1.0 mJy bcam~^. The lo- 
cal noise near the brightest three sources is reduced to 
< 2.0 mJy beam~^. 

For an observation of ^ 10 hours, the theoretical ther- 
mal noise level'^ for the GMRT at 153 MHz is estimated 
to be ^ 0.2 — 0.3 mJy beam~^. For our observation, 
the measured noise of ^ 1.0 mJy beam~^ in the cen- 
tral part of the field is a factor of 3 to 5 larger. In ap- 
parently source-free parts of the outer imaged area this 
drops to 0.75 mJy beam"-'^, indicating that residual side- 
lobe structure has a noticable contribution to the central 
noise. Although significantly improved as compared to the 
self-calibrated case, the residual sidelobe structure near the 
three bright sources, all located near the half-power ra- 
dius of the primary beam area, limits the local dynamic 
range to ^ 600. We expect that the residual sidelobes are 
caused by direction-dependent visibility amplitude errors 
due to pointing errors and non-circular primary beam pat- 
terns (e.g., Bhatnagar et al. 2008). 

Another suspect for adding to the overall noise is resid- 
ual RFI (see also Section 3.3). Especially on the most af- 
fected, short baselines, it can be difficult to distinguish be- 
tween persistent broad-band RFI and true sky signal. In our 
approach we have attempted to save as many short base- 
lines as possible, to preserve the sensitivity to large-scale 
emission. This could explain why our central noise level 
of ~ 1.0 mJy beam~^ is slightly higher than the quoted 
0.7 mJy beam~^ found by Sirothia et al. (2009); Ishwara- 
Chandra et al. (2010) for similar deep fields. In general, 
we find that a straightforward comparison of noise levels 
is non-trivial, as it depends on the flagging & data selec- 
tion strategy, the visibility weighting scheme used during 
imaging, and on how and where in the image the noise is 
measured. 

3. Catalog Construction 

In this section we describe the construction of the source 
catalog that was extracted from the 153 MHz image of the 
Bootes field. The image after SPAM calibration was scaled 
down by 4 percent to incorporate a correction to the flux 
calibrator scale (sec Section 3.3). Next, the image was cor- 
rected for primary beam attenuation with a circular beam 
model"* out to a 1.9 degree radius (out to 0.3 times the cen- 
tral beam response, slightly beyond the half-power radius) , 
yielding a survey area of 11.3 square degrees. Figure 1 shows 
a map of the local noise of the circular survey area, which 

^ Derived using the GMRT User Manual; http://gmrt.ncra.tifr. 
res. i n/gm rt_h page/Use rs/doc/manua I/Use rsManua I/index, htm I 

See the GMRT User Manual 
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BACKGROUND RMS [mjy/beam] 



12 3 4 




1440 38 36 34 32 30 28 26 24 

RIGHT ASCENSION (J2000) 



Figure 1. The grayscale map represents the local noise 
as measured in the (primary beam corrected) Bootes 
field image. The overplotted contours mark lines of equal 
RMS at [1.0,1.4,2.0,2.8,4.0] mjy beam"^ (white above 
2.0 mJy beam~^). The local enhancements in RMS coincide with 
the positions of the brightest sources. 



has a central area with a noise level of ~ 1.0 mJy beam~^, 
a global increase of the local noise to 2.0-2.5 mJy beam^^ 
when moving towards the field edge, and several small ar- 
eas around bright sources where the local noise is approxi- 
mately twice the surrounding noise level. The overall noise 
was found to be 1.7 mJy beam~^. 

Although we used multi-frequency synthesis, the 
0.25 MHz width of individual frequency channels causes 
bandwidth smearing during imaging. Applying the stan- 
dard formula (e.g., Thompson 1999) to our case, the radial 
broadening of sources at the edge of the field (at 1.9 degree 
from the pointing center) is estimated to be ~ 11", which 
is half the minor axis of the restoring beam. Similarly, the 
visibility time resolution of 16.8 sec causes time-average 
smearing in the order of ~ 8" at the field's edge. This may 
appear problematic, but the majority of sources are de- 
tected within the inner part of the primary beam (both 
effects scale linearly with radial distance from the pointing 
center). The total flux density of smeared sources is con- 
served, but the peak fiux drops. The effect of smearing on 
source extraction is discussed in Section 3.2. 

3.1. Source Extraction 

We used the BDSM software package (Mohan 2008) to ex- 
tract sources from our image. With our settings, BDSM 
estimates the local background noise level (7^ over the map 
area, searches for pixels > 5 CTl , expands the 5 detections 
into islands by searching for adjacent pixels > Sctl, reject- 
ing islands smaller than 4 pixels. BDSM fits the emission in 



each island with one or multiple gaussians, combines signif- 
icantly overlapping gaussians into sources, and determines 
the fiux densities, shapes and positions of sources (includ- 
ing error estimates, following Condon 1997). This approach 
has led to very few false source detections (see Section 3.2). 

BDSM detected 644 islands, for which the 935 fitted 
gaussians were grouped into 696 distinct sources. Of these, 
499 sources were fitted with a single gaussian. Visual in- 
spection of the image, complemented with a comparison 
against a very deep 1.4 GHz map (see Section 3.3), re- 
sulted in the removal of 16 source detections in the near 
vicinity of six of the seven brightest sources. We also re- 
moved 4 sources that extended beyond the edge of the im- 
age. To facilitate total flux density measurements at the 
high fiux end, we combined multiple source detections in 
single islands, and manually combined 50 additional source 
detections that were assigned to different islands but ap- 
peared to be associated in either the 153 MHz image or the 
deep 1.4 GHz map. The combined flux density is the sum 
of the individual components, while the combined position 
is a flux- weighted average (centroid). Error estimates of the 
positions and flux densities of the components are propa- 
gated into error estimates of the combined flux density and 
position. The final catalog consists of 598 sources. 

Table 3 contains a sample of the full source catalog avail- 
able in the CDS^ database. The source entries include the 
corrections discussed in the remainder of Section 3. The 
columns are: (1) source identifier, (2,3) centroid position 
right ascension and declination, (4,5) centroid position un- 
certainty along the right ascension and declination direc- 
tions, (6,7) integrated source fiux and uncertainty, (8) local 
background RMS noise level, and (9) number of gaussians 
fitted to the source. 

3.2. Completeness and Contamination 

We estimated the completeness of the 153 MHz catalog by 
performing Monte-Carlo simulations. The source extraction 
process generated a residual image from which all detected 
source flux was subtracted. For our simulation, we inserted 
1000 artificial point sources into the residual image, and 
used the same mechanism as described in Section 3.1 to 
extract them. The artificial source positions were selected 
randomly, but never within 50" of another source, a blanked 
region (near the image edge) or a > 10 mJy residual. The 
source peak fluxes S were chosen randomly within the range 
3 mJy to 3 Jy, while obeying the source count statistic 
dN/dS oc S~^'^, which produced statistically sufficient de- 
tection counts in all logarithmic fiux bins. The simulation 
was repeated 20 times to improve the accuracy and to de- 
rive error estimates. The detection fraction as a function of 
peak flux is plotted in Figure 2. 

In this approach we have ignored several effects that 
may influence the detectability of sources, such as (i) the 
intrinsic size of sources, (ii) bandwidth- and time-average 
smearing, (iii) calibration errors, and (iv) imaging and de- 
convolution. The source detection algorithm uses a peak de- 
tection threshold. Sources that are resolved, smeared due to 
bandwidth- and time-averaging, or defocussed due to cali- 
bration errors will therefore have a decreased probability of 
detection. We have not attempted to model for the angu- 
lar size distribution of sources at this frequency. However, 

http://cdsweb.u-strasbg.fr/ 
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ID 


RA 


DEC 


ARA 


ADEC 


S153 


A5i53 


Noise 


gauss 




[1 


[1 


["] 


n 


[mJy] 


[mJy] 


[mJy beam^^] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


J142632+350815 


216 63367 


35 13754 


5.2 


4.3 


649.7 


65.1 


2.8 


1 


J143810+340459 


219 54576 


34 08315 


3.9 


2.8 


641.2 


64.2 


2.0 


1 


1143850-1-330225 


219 71068 


33 04036 


10.0 


8.8 


577.8 


58.2 


2.9 


2 


1143134-1-351505 


217 89466 


35 25165 


5.2 


3.5 


567.8 


56.9 


1.6 


2 








71 1 


fiO 


54"^ n 




9 5 


7 


J142445+341816 


216.19020 


34.30447 


7.2 


5.3 


507.1 


50.9 


2.9 


1 


J1430554-341933 


217.73287 


34.32597 


3.9 


3.0 


505.0 


50.6 


1.2 


2 


J143256-h353339 


218.23353 


35.56104 


4.7 


4.0 


503.8 


50.5 


2.0 


1 


J143930-h344628 


219.87580 


34.77470 


9.2 


7.3 


489.8 


49.4 


2.5 


2 


J143128+355252 


217.86699 


35.88132 


11.6 


11.6 


459.4 


46.5 


2.3 


3 



Table 3. Sample of the full source catalog of GMRT 153 MHz sources. The columns are described in the text. 
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Figure 2. Detection fractions and completeness as a function of flux. Left: Result from 20 Monte-Carlo simulations, in which 
1000 point sources with varying fluxes were inserted into the residual image, followed by a source extraction. The horizontal axis 
denotes the input (peak) flux. The solid line is the average detection fraction, and the dotted lines denote the 1 a uncertainty. 
Middle: Result from scaling the peak flux and total flux density of a subset of high S/N sources with varying sizes down to the 
detection limit, and applying the 5o"l detection criteria. The dashed line is the detection fraction as a function of peak flux, the 
solid line is the detection fraction as a function of total flux density, and the dotted line is the dashed line shifted in flux by 
15 percent to approximately match the solid line. Right: Estimated completeness of the extracted source catalog as a function of 
(total) flux limit. 



previous observations show that the major fraction of low- 
frequency sources are unresolved at ~ 25" resolution (e.g., 
Cohen et al. 2004; George & Ishwara-Chandra 2009). In our 
catalog, more than 90 percent of the sources appear to have 
simple, near-gaussian morphologies. Assuming the angular 
size distribution of sources changes slowly with source flux 
density, and assuming that smearing and defocussing af- 
fects sources of varying brightness in a statistically equal 
way, we can estimate the resolution bias from the catalog 
itself. For this purpose, we select a subset of 214 high S/N 
sources with peak fluxes between 12 and 20 Cl and sim- 
ple morphologies. The flux densities of these sources were 
scaled down by a factor of 4 to create an artificial popula- 
tion of sources near the detection threshold. After applying 
the 5 detection criterium, we determined the detection 
fractions, both as a function of peak flux and total (in- 
tegrated) flux density (Figure 2). Although this approach 
suffers from low number statistics, the general trend of both 
detection fraction functions is similar but shifted upwards 
in flux by ~ 15 percent. We therefore assume that the total 
flux detection fraction is approximated by the peak flux de- 



tection fraction derived from the Monte-Carlo simulations, 
shifted upwards in flux by 15 percent. 

To estimate the completeness of the catalog at various 
flux limits, the total flux detection fraction must be mul- 
tiplied by the (normalized) true source flux distribution, 
and integrated from the flux limit upwards. The Euclidean- 
normalized differential source counts in Section 4.1 approx- 
imately follow a power-law over the larger part of the flux 
range of our sample, with a slope of 0.91. The source flux 
distribution over the same flux range can therefore be ap- 
proximated by a power-law with slope of 0.91 — 5/2 = 
— 1.59. As it involves very few counts, the deviation of the 
power-law from the true source counts at the high flux 
end has little effect on our estimates, and is therefore ig- 
nored. Figure 2 shows the resulting completeness estimate 
as a function of total flux. From this plot, the estimated 
completeness is 70 percent at 5.1 mJy, ^ 92 percent at 
10 mJy and > 99 percent at 20 mJy. 

A known effect that arises from the dcconvolution pro- 
cess is CLEAN bias (e.g., Condon et al. 1994, 1998; Becker 
et al. 1995). This is a systematic negative offset in the 
recovered flux densities after dcconvolution, probably the 
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result of false CLEANing of sidelobe peaks in the dirty 
beam pattern. One can estimate the CLEAN bias by in- 
jecting artificial sources into the visibility data and compare 
the recovered flux densities after imaging & deconvolution 
with the injected flux densities. We have not attempted 
this approach, but instead taken precautions to minimize 
the CLEAN bias effect. In our case, the dirty beam is 
well-behaved due to a relatively uniform UV-coverage from 
two extended observing runs, in combination with multi- 
frequency synthesis and a robust weighting parameter of 
0.5 (slightly towards natural weighting). CLEAN bias is 
further suppressed through the use of tight CLEAN boxes 
in the imaging & deconvolution process. 

For an estimate of the contamination of the catalog with 
fake detections, we compare the GMRT 153 MHz image and 
extracted source catalog against the results from the deep 
WSRT 1.4 GHz survey of the Bootes field by de Vries et al. 
(2002). The 153 MHz and 1.4 GHz observations are weh 
matched in terms of survey area and resolution 7 square 
degrees and 13" x 21" for WSRT, respectively). The typ- 
ical noise over the WSRT survey area is 28 fiJy beam~^. 
For a spectral index of —0.8, the 1.4 GHz observations are 
^ 10 times more sensitive. We restrict our comparison to 
a 1.4 degree radius circular area to avoid the noisy edge 
of the deep 1.4 GHz survey. For all of the 399 sources de- 
tected at 153 MHz we find a counterpart in the 1.4 GHz 
map (383 sources were automatically matched within a 25" 
search radius, while the remaining fraction of sources with 
complex morphology were confirmed manually). We could 
not match the full GMRT area, but considering that our 
source extraction is based on local noise and that the false 
detections only appeared to occur near a few bright sources, 
we estimate that the contamination of our complete catalog 
over the full survey area is < 1 percent. 

3.3. Astrometric and Flux Uncertainty 

For an estimate of the astrometric uncertainty, we compare 
the source positions in the GMRT 153 MHz catalog against 
catalog source positions from the deep WSRT 1.4 GHz map 
of de Vries et al. (2002). For our position comparison, we 
only use sources whose flux profile is accurately described 
by a single gaussian, and whose peak flux Sp is at least 
10 (Tl. This bypasses most of the position errors that arise 
from low signal-to-noise (S/N, or 5^/(7^), different grouping 
of gaussians and spectral variations across sources. Using 
a search radius of 10", we cross-match 126 sources in both 
catalogs. From this, we measure a small mean position off- 
set in RA and DEC of (Aa, AJ) = (0.11", 0.09"). We cor- 
rect the catalog positions for this small offset. The esti- 
mated RMS scatter around this offset is % ^ = 1.32". The 
S/N- independent part of the positional uncertainty of the 
1.4 GHz sources is 0.44" (de Vries et al. 2002), therefore we 
derive an absolute astrometric uncertainty for the 153 MHz 
sources of 1.24". We quadratically add this to the calculated 
(S/N-dependent) position accuracies in the 153 MHz source 
catalog (see Section 3.1). 

The uncertainty of the flux scale transferred from the 
calibrator 3C 286 to the target Bootes held is influenced 
by several factors: (i) the quality of the calibrator obser- 
vational data, (ii) the accuracy of calibrator source model, 
and (iii) the difference in observing conditions between the 
calibrator and target field. Because of the relatively large 
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Figure 3. Model of the radio spectrum of 3C 286. The mea- 
surements (dots) above 153 MHz (dashed line) were used to fit 
the original Perley-Taylor 1999.2 model (dotted line), which ap- 
pears to overestimate the flux density at 153 MHz. By adding 
one measurement at 74 MHz, we fit an alternative model (solid 
line) that may give a better prediction at 153 MHz. 



uncertainty in the flux scale at low frequencies, we discuss 
in some detail the issues that influence these factors. 

The quality of the calibrator data is most noticeably 
affected by RFI and by ionospheric phase rotations. The 
repeated observation of 3C 286 every ^ 45 minutes during 
the observing enabled us to monitor these effects over time. 
The mild fluctuations in the initial (short interval) calibra- 
tion gain phases at the start of the data reduction showed 
that the ionosphere was very calm during both observ- 
ing nights, therefore we exclude the possibility of diffrac- 
tion or focussing effects (e.g., Jacobson & Erickson 1992). 
Apparent flux loss due to ionospheric phase rotations was 
prevented by applying the (short interval) gain phase cor- 
rections before bandpass- and amplitude calibration (see 
Section 2.2). 

RFI was continuously present during both observing 
sessions. This mainly consisted of persistent RFI over the 
full band, most noticeably on the shortest (central square 
and neighbouring arm antenna) baselines, and of more 
sporadic events on longer baselines during one or more 
time stamps and/or narrow frequency ranges. The sporadic 
events were relatively easy to recognize and excise, but for 
the persistent RFI this is much more difficult due to a lack 
of contrast between healthy and affected data on a single 
baseline. Some of the shortest, most affected baselines were 
removed completely. On longer baselines, persistent RFI 
from quasi-stationary sources can average out due to fringe 
tracking (Athreya 2009), but does add noise. Large mag- 
nitude RFI amplitude errors in the visibilities may result 
in a suppression of the gain amplitude corrections. Because 
these effects are hard to quantify, we adopt an ad-hoc 2 per- 
cent amplitude error due to RFI. 

Because 3C 286 is unresolved (< 2.5") within a 20" to 
25" beam, we used a point source for the calibrator model 
with a flux density of 31.01 Jy (Section 2.2). The utilized 
Perley-Taylor flux density at 153 MHz is an extrapolata- 
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Frequency 
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VLSS'' 


74 


30.26 ±3.08 


34.67 


30.26 


dC 


151 


26.31 


31.25 


on o 1 

29.81 




151 


26.53 


31.25 


29.81 




159 


30.0 ±7.0 


30.94 


29.65 




178 


24.0 


30.22 


29.26 


WENSS'' 


325 


27.12 ± 1.63 


25.96 


26.11 



" Perley- Taylor 

Cohen et al. (2007) 
"Hales et al. (1988) 

Hales et al. (2007) and references therein. Note that we use the peak flux 

rather than the integrated flux density, as 3C 286 is unresolved on the 7C 

resolution. 
" Edge et al. (1959); Bennett (1962) 
^ Pilkington & Scott (1965); Gower et al. (1967) 

* Rengelink et al. (1997). Note that we use the peak flux rather than the 
integrated flux, as 3C 286 is unresolved on the WENSS resolution. 

Table 4. Flux measurements of 3C 286 from different radio survey catalogs at low 
frequencies. 



tion of VLA flux density measurements at 330 MHz and 
higher, using a fourth order polynomial in log- log space. 
In Table 4 a comparison is presented between flux density 
measurements of 3C 286 in various sky surveys at low fre- 
quencies and the predicted flux densities from the Perley- 
Taylor model. Although there is a large variation in the flux 
differences at the different frequencies, there appears to be a 
overestimation by the Perley-Taylor model below 200 MHz 
due to a spectral turnover of 3C 286 below ~ 300 MHz. For 
this reason, we re-fitted the polynomial to the original data 
points plus the additional 74 MHz VLSS measurement, as- 
suming 5 percent errors for all data points (Figure 3). Our 
new model is given by: 

logio(^J = 1-24922 - 0.434710 logio(i^) 

- 0.174786(logio(j/))2 + 0.0251542(loglo(^/))^ (1) 

with V the frequency in GHz and S^, the flux density in Jy. 
This parametrization results in a slightly better fit to the 
flux densities from the various surveys (Table 4), although 
the large scatter remains. For the center of the GMRT band 
at 153.1 MHz, the new model predicts a flux density of 
29.77 Jy. We have adopted this new model by scaling the 
image by the ratio 29.77/31.01 before primary beam cor- 
rection (see start of Section 3). The accuracy of this flux 
scale is closely related to the accuracy of the 74 MHz data 
point, for which a 10 percent flux error is given. Assuming 
the error in the 330 MHz is much smaller, and 153 MHz is 
roughly half way from 74 to 330 MHz in logarithmic fre- 
quency, we set an upper limit of 5 percent error on the 
adopted flux scale of 3C 286 at 153 MHz. 

The presence of other sources in the 3.1 degree FoV 
around 3C 286 further complicates the flux calibration. 
For example, the total apparent flux density in the Bootes 
field that was extracted through CLEAN deconvolution is 
^ 46 Jy, which may be typical lower limit for any blind 
field. If this flux is distributed over many sources that are 
individually much fainter than the calibrator, then the net 
effect of this additional flux is only noticeable on a small 
subset of the shortest baselines, while calibration utilizes all 
basehnes. Inspection of the 3C 286 field at 74 MHz (VLSS; 
Cohen et al. 2007) and 325 MHz (WENSS; Rengelink et al. 



1997) does identify two relatively bright sources within 
0.7 degrees of 3C 286 with estimated apparent GMRT 
153 MHz flux densities of 5.3 and 2.7 Jy, respectively. These 
sources will cause a modulation of the visibility amplitudes 
across the UV-plane. We performed a simple simulation, 
in which we replaced the measured 3C 286 visibilities with 
noise-less model visibilities of three point sources, being 
3C 286 and the two nearby sources, and calibrated these 
visibilities against a single point source model of 3C 286, 
using the same settings as in the original data reduction. 
We found that the combined gain amplitude for all anten- 
nas and all time intervals was 1.000 ±0.004. For individual 
10 minute time blocks, the largest deviation from one was 
1.00 percent, which indicates the magnitude of the possible 
error when using a single 10 minute calibrator observation 
on 3C 286. The small deviations per time block are trans- 
ferred to the target field data, but these are suppressed 
by amplitude self-calibration against the target field source 
model. We set an upper limit of 1 percent due to the pres- 
ence of other sources in the FoV of the calibrator. 

While transferring the flux scale from calibrator to tar- 
get field, the derived gain amplitudes need to be corrected 
for differences in sky temperature due to galactic diffuse ra- 
dio emission (e.g., Tasse et al. 2007), which is detected by 
individual array antennas but not by the interferometer. 
The GMRT does not implement a sky temperature mea- 
surement, therefore we need to rely on an external source of 
information. From the Haslam et al. (1982) all-sky map, we 
find that the mean off-source sky temperatures at 408 MHz 
as measured in the 3C 286 and Bootes field are both ap- 
proximately 20 ± 1 degree. Applying the formulae given by 
Tasse et al. (2007) for the GMRT at 153 MHz^, we estimate 
a gain inaccuracy of ~ 2 percent at most. 

Another effect that may influence the gain amplitude 
transfer between calibrator and target field is an elevation- 
dependent gain error. This is a combination of effects such 
as structural deformation of the antenna, atmospheric re- 
fraction and changes in system temperature from ground 

® We adopted the GMRT system parameters from 
http://www.gmrt.ncra.tifr.res.in/gmrt_hpage/Users/doc/manual/ 
UserslVlanual/nodel3.html 
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radiation. According to Chandra ct al. (2004), the effect 
on amplitude is rather small, if not negligible, for GMRT 
frequencies of 610 MHz and below. Furthermore, the rela- 
tively short angular distance of 13.4 degrees between 3C 286 
and the target field center causes the differential elevation 
error to be limited. Elevation dependent phase errors are 
not relevant, because we don't rely on the calibrator to re- 
store the astrometry. For our observation, we assume that 
elevation-dependent effects can be ignored. 

To incorporate the effects discussed above plus some 
margin, we quadratically add a relative flux error of 10 per- 
cent to the flux density measurement errors in the source 
catalog. 



4. Analysis 

The 598 radio sources in our 153 MHz catalog form a statis- 
tically significant set, ranging in flux density from 5.1 mJy 
to 3.9 Jy. Because of the large survey area, cosmic variance 
is expected to be small. Radio images at 153 MHz of a se- 
lection of extended sources are presented in Appendix A. 
In this section, we discuss two characteristics of the sur- 
vey: number counts and spectral indices. Further analyses 
will be done in a subsequent article. When properly cor- 
rected for incompleteness, the number counts are an objec- 
tive measure of the 153 MHz source population, that can be 
compared against models and other surveys. For individual 
sources, the spectral index can help to classify sources, and 
identify rare sources such as HzRGs and cluster halos & 
relics. 

4.1. Differential Source Counts 

We derived the Euclidean-normalized differential source 
counts from the catalog. Because the source extraction cri- 
teria vary over the survey area, we used Figure 2 to cor- 
rect for the missed fraction per flux bin. This mainly af- 
fects the lowest two bins. Furthermore, the combined ef- 
fect of random peaks in the background noise and a peak 
detection criteria causes a selection bias for positively en- 
hanced weak sources (Eddington bias). In general, noise 
can scatter sources into other flux bins, most noticeably 
near the detection limit. Our attempts to correct for this 
effect through Monte-Carlo simulations were numerically 
unstable due to the low number counts in the lowest flux 
bin. The effect on the higher bins was minimal, therefore 
we omit the Eddington bias corrections. The Euclidean nor- 
malized differential source counts are tabulated in Table 5. 
Figure 4 shows our source counts in comparison with the 
(little) other observational data that is available for this fre- 
quency at comparable sensitivity, as well as source counts 
derived from two models. 

At the high flux end (0.1 — 1 Jy), there is a good agree- 
ment between our source counts and those derived from 
part of the 7C survey at 151 MHz (McGilchrist et al. 
1990). Ishwara-Chandra et al. (2010) present source counts 
of 765 sources in a single-field, deep GMRT 153 MHz obser- 
vation similar to ours, but with over twice the bandwidth 
(16 MHz). Their source counts go down to 2.5 mJy. While 
their source counts are roughly equal to ours at the high 
flux end, they are increasingly lower towards lower flux lev- 
els. Their source counts are well approximated by a single 
power- law slope of 1.01, which is steeper than the value of 




flux [Jy] 



Figure 4. Euclidean-normalized differential source counts for 
tlie full 153 MHz catalog of 598 sources, distributed over 14 log- 
arithmic flux bins ranging from 4.75 mJy to 3 Jy (black dots; we 
have omitted the highest flux bin, containing just one count). 
Also plotted are 151 MHz source counts by McGilchrist et al. 
(1990), part of the 7C catalog (black triangles), 153 MHz source 
counts by Ishwara-Chandra et al. (2010, white triangles), model 
151 MHz source counts by Jackson (2005) for a ACDM cosmol- 
ogy (dotted line) and an Qm = 1 cosmology (dash-dotted line), 
a 151 MHz source count model by Wilman et al. (2008, solid 
line), and the power-law slope of 0.91 (dashed line) that fits our 
data points best. 
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0.0945 
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0.1504 


36 
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0.2383 


19 


1375.27 ±315.51 


0.3777 


16 
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0.5986 


11 


3169.75 ±955.72 


0.9487 
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3449.72 ± 1408.34 


1.5036 


7 


8030.29 ±3035.16 


2.3830 
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2288.93 ±2288.93 



Table 5. Euclidean-normalized differential source counts 
(including poissonian error estimates) for the full 153 MHz 
catalog of 598 sources, distributed over 14 logarithmic flux 
bins ranging from 4.75 mJy to 3 Jy. 



0.91 from our data over the range 40-400 mJy. This sug- 
gests a non-linear flux scale difference. We can offer little 
explanation on this apparent discrepancy. Within 2-3 times 
the (Icr) error bars, the two source count results may still 
be considered consistent. 

George & Stevens (2008) have determined source counts 
from a shallower GMRT 153 MHz field (central RMS 
3.1 mJy beam""'^) centered around e Eridanus. From bin- 
ning 113 source fluxes over a range of 20 mJy-2 Jy they 
derive a single power-law slope of 0.72 (not plotted). The 
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most likely explanation for this much flatter slope is the 
large uncertainties in their counts due to low-number statis- 
tics. 

Jackson (2005) constructed a 151 MHz source count 
model, based on an extrapolation of source counts > 0.2 Jy 
from the 3CRR and 6C catalogs (Laing et al. 1983; Hales 
et al. 1988). Two cosmological scenarios are considered, 
namely an ft^ — 1 cosmology and a ACDM cosmology, 
the latter being today's generally accepted cosmology (e.g., 
Komatsu et al. 2011). In this model, FRii radio sources 
dominate the counts above ~ 50 and ~ 20 mJy, for the two 
cosmologies respectively. Below, FRi sources are the most 
dominant population down to below our detection thresh- 
old, which causes a flattening of the counts below ~ 20 mJy 
in both cosmologies. 

Figure 4 shows that the source count models for both 
scenarios roughly match our observed source counts near 
500 mJy, but increasingly underestimate the counts towards 
lower flux densities. Between 20 and 200 mJy, the approxi- 
mately constant model slope for both scenarios is 0.98 and 
1.19, respectively, steeper than the value of 0.91 derived 
from our data. The observed source counts shows no clear 
evidence of flattening towards lower flux levels. 

Wilman et al. (2008) have generated model source 
counts at 151 MHz from a semi-empirical (ACDM) cos- 
mological simulation, using (extrapolated) radio luminos- 
ity functions for several different populations of sources. In 
this model, FRii radio sources dominate the counts above 
~ 200 mJy, while below FRi sources dominate down to 
beyond our detection threshold. Figure 4 shows an ap- 
proximate agreement between the model and our obser- 
vations, although the average power-law slope between 10 
and 200 mJy of 0.83 is flatter than the value of 0.91 derived 
from our data. 

For the models that are discussed here we find that all 
predicted source counts deviate from the observed source 
counts. The model source count slopes appear to be influ- 
enced mostly by the assumed contribution of FRi sources. 
On the other hand, the observational results from Ishwara- 
Chandra et al. (2010) also deviate from our result. More 
source counts from similar deep surveys at this frequency 
are needed to put stronger constraints on the source count 
models. Preliminary source counts from a deep GMRT 
153 MHz survey three times the area of our current survey 
(> 1000 sources; Intema et al., in preparation) match up 
closely with our current results. And future source counts 
from the ongoing GMRT 153 MHz sky survey (TGSS) 
should provide a more robust reference to compare the 
source count models and our results against. 

4.2. Spectral Indices 

Because of the good match in resolution between the 
GMRT 153 MHz image and the deep WSRT 1.4 GHz im- 
age from de Vries et al. (2002), we can accurately deter- 
mine spectral indices over a decade in frequency. The survey 
depths at 153 MHz and 1.4 GHz are equal for sources with 
a spectral index of —1.6. Due to the high detection rate 
of 1.4 GHz sources at 153 MHz positions (Section 3.2), we 
do an automated search for 1.4 GHz counterparts within 
25" of the 153 MHz sources and ignore all sources for 
which we don't find counterparts. The spectral indices of 
417 matched sources arc plotted in Figure 5. We find a 
median spectral index of —0.76, which is similar to the me- 
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Figures. Spectral index between 1.4 GHz and 153 MHz for 
417 sources, plotted as a function of the 153 MHz flux (black 
dots). The detection limit is almost fully determined by the 
153 MHz survey due to the sensitivity of the 1.4 GHz ob- 
servations. The median spectral index is —0.76 (dashed line). 
16 sources have a spectral index below —1.3 (dotted line). The 
spectral indices are binned to emphasize the spectral flattening 
towards lower flux densities (white squares; see Table 6). 



dian values of -0.79 (Cohen et al. 2004), -0.85 (Ishwara- 
Chandra & Marathe 2007), -0.82 (Sirothia et al. 2009) and 
—0.78 (Ishwara-Chandra et al. 2010) found for similar high- 
resolution, low- frequency surveys at 74 and 153 MHz. The 
small differences are most likely caused by differences in the 
completeness limits of the catalogs, as the median spectral 
index varies with flux density (see below). 

Figure 5 also plots the median spectral index in 6 log- 
arithmic flux bins, which clearly highlights the flattening 
trend of the mean spectral index towards lower flux densi- 
ties. The unique combination of our deep 153 MHz catalog 
and the very deep 1.4 GHz catalog makes that the median 
spectral index is unbiased down to the lowest 153 MHz flux 
densities. The median spectral index is ^ —0.9 for > 0.5 Jy 
sources, and flattens to ~ —0.7 for < 50 mJy sources. A 
similar trend appears to be present in the spectral index 
distribution by Ishwara-Chandra et al. (2010) for sources 
between 153 MHz and various frequencies (610 MHz and 
higher), but may be biased by catalog flux limits at the 
higher frequencies. A similar flattening trend is also seen 
at both lower and higher frequencies, e.g., between 74 MHz 
and 1.4 GHz by Cohen et al. (2004) and Tasse et al. (2006), 
between 1.4 GHz and 325 MHz by de Vries et al. (2002), 
between 325 MHz and 1.4 GHz by Zhang et al. (2003) and 
Owen et al. (2009), between 610 MHz and 1.4 GHz by Bondi 
et al. (2007), and between 1.4 and 5 GHz by Prandoni et al. 
(2006). 

Recently, Bornancini et al. (2010) found that, for a near- 
complete sample of radio galaxies, there is no evidence for 
spectral steepening or flattening due to redshifted curved 
radio spectra. In fact, most of their sources have straight 
power-law spectra from 74 MHz to 4.8 GHz. This sug- 
gests that the observed spectral flattening towards lower 
flux densities results dominantly from a correlation between 
source luminosity and spectral index {P — a correlation). 
This correlation is found to exist for FRii galaxies (e.g.. 
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Flux bin center 

[Jy] 


Counts 


Spectral index" 


0.0083 


133 


-0.710 ±0.284 


0.0232 


133 


-0.709 ±0.208 


0.0643 


84 


-0.784 ±0.188 


0.1786 


38 


-0.847 ±0.173 


0.4966 


19 


-0.897 ±0.136 


1.3799 


10 


-0.908 ±0.123 



" Median ± RMS around median. 

Table 6. Median spectral index between 1.4 GHz 
and 153 MHz for 417 sources in 6 logarithmic 
(153 MHz) flux bins ranging from 5.0 mjy to 
2.3 Jy. 



Blundell et al. 1999). According to the models by Jackson 
(2005) and Wilman et al. (2008), FRii galaxies are the dom- 
inant 153 MHz source population in our sample at higher 
flux levels (> 20 — 200 mJy), which can explain part of the 
flattening observation. For lower flux levels, the flattening 
trend continues, which suggests that either FRii galaxies 
are still a significant source population, or the FRi galaxies 
dominate and also follow a, P — a correlation. 

From Figure 5 we highlight a small group of 16 USS 
sources that have a spectral indices lower than —1.3. The 
lowest spectral index is —1.65, which indicates that the 
sample selection limit is indeed dominated by the 153 MHz 
survey depth. Table 7 lists these sources in decreasing flux 
order. Despite the relatively large uncertainty in 153 MHz 
flux density for the faintest sources, the spectral index is 
still relatively well constrained due to the large frequency 
span. 

The fraction of USS sources is 3.8 percent (16 out of 
417), which is mostly dependent on the survey detection 
limit at 153 MHz. It is not straightforward to compare this 
fraction with the few other 153 MHz surveys, as they use 
different high frequencies, with different detection limits at 
the high- and low-frequency end. Using the same criteria, 
Sirothia et al. (2009) find a USS fraction of 3.7 percent (14 
out of 374). Applying the same criteria to the source catalog 
by Ishwara-Chandra et al. (2010) yields a USS fraction of 
3.0 percent (19 out of 639). These three fractions appear to 
have some consistency. 

The angular distribution of the USS sources is quite 
peculiar. There are 6 sources that form 3 pairs within 6' 
of each other. These pairs, together with the remaining 
10 single sources, appear to be randomly distributed across 
the FoV. Visual inspection of the pairs in the images at 
153 MHz and 1.4 GHz did not reveal any obvious artifacts 
in the background near these sources, which makes it less 
likely that these sources are fake detections. There is no 
visual evidence in the radio maps that the pair compo- 
nents are physically connected. Through Monte-Carlo sim- 
ulations we determined that, out of 16 sources with ran- 
dom positions over a 1.5 degree radius field (which is ap- 
proximately the 1.4 GHz field radius), the chance of finding 
3 pairs within 6' is 1.4 ±0.4 percent. Therefore, it is unlikely 
that these pairs appeared by chance. Further investigation 
is needed to establish the true nature of these pairs. 

Croft et al. (2008) examine 4 candidate HzRGs (which 
they labelled A, B, C and E) based on their steep (< —0.87) 
spectral index between 1.4 GHz and 325 MHz. As their can- 
didate sources are also present in the 153 MHz catalog, we 



complement their data with our new spectral index mea- 
surements in Table 8. Both sources A and E appear to have 
fairly straight power-law spectra down to 153 MHz, while 
sources B and C appear to undergo considerable spectral 
flattening. Based on our selection criteria for USS sources, 
only source A would be considered a candidate HzRG. 



4.3. Identification Fraction of Radio Sources versus Spectral 
Index 

We investigate the NIR A'-band identification fraction of 
radio sources as function of spectral index (e.g., Wieringa 
1991). This links together two known correlations, namely 
the correlation between spectral index and redshift (Tielens 
et al. 1979; Blumcnthal & Milcy 1979) and the correlation 
between A'-band magnitude and redshift {K—z correlation; 
e.g., Willott et al. 2003; Rocca-Volmerange et al. 2004). 
We remark that the correlation between spectral index and 
redshift is not a tight correlation; for a given spectral index 
there is a range of possible redshifts for a radio source (e.g., 
Miley & De Breuck 2008). But on average, steep spectrum 
sources are located at higher redshifts, and are consequently 
(on average) more difflcult to detect in optical and NIR 
bands. Using A'-band has an advantage over (also available) 
shorter wavelength bands as this band suffers the least from 
extinction. 

We have identified possible optical counterparts of 
the radio sources using the FLAMEX AT^-band catalogue 
(Elston et al. 2006). This survey covers 7.1 square degrees 
within the Bootes field. For the NIR identification we use 
the likelihood ratio technique described by Sutherland & 
Saunders (1992). This allows us to obtain an association 
probability for each NIR counterpart, taking into account 
the NIR magnitudes of the possible counterparts. Before 
the cross identification, we removed all radio sources lo- 
cated outside the coverage area of the A'^-band images, or 
located within 20 pixels of the edge or other blanked pixels 
within individual A'^-band frames. 

Following Sutherland & Saunders (1992) and Tasse 
et al. (2008), the probability that a NIR counterpart with 
magnitude m is the true NIR counterpart of the radio 
source is given by the likelihood ratio 



LR(r, m) 



0{< m)cxp {-r'^/2) 
27r (Tq ffg p{< m) 



(2) 



with m the AT^-band magnitude of the NIR candidate, 
0{< m) the a priori probability that the radio source has 
a NIR counterpart with a magnitude brighter than m, and 
p{< m) the surface number density of NIR sources with a 
magnitude smaller than m. (t| and ct^ are the quadratic 
sums of the uncertainties in the radio and NIR positions in 
right ascension and declination, respectively. For the radio 
source positions and uncertainties we have taken the val- 
ues from the 1.4 GHz WSRT catalogue, as the astrometric 
precision is better than the 153 MHz GMRT data. For the 
FLAMEX survey we have adopted 0.3" for the position 
uncertainty (this is half the CCD pixel size, as source posi- 
tions are truncated to pixel positions). The uncertainty- 
normalized angular distance is given by the parameter 
r ^ -y(A„/cr„)2 ^ (A^-/o-^)2, with A the positional differ- 
ence in either a oi 6 between the possible NIR counterpart 
and radio source. 
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ID'' RA*" DEC*" 153 MHz flux 1.4 GHz flux 



[mJy] [mJy] 



J142656+352230 


14i^26'"56, 


.46" 


35°22'30, 


.8" 


264, 


,2 ±26.6 


11.31 ±0.46 


-1 


.43 


± 


0, 


,05 


J143506+350059 


14i^35'"06, 


.89" 


35°00'58, 


.2" 


142, 


,8 ± 14.5 


4, 


.02 ± 





.16 


-1 


.62 


± 


0, 


,05 


J143118+351549 


14'^31'"18 


.29" 


35°15'49. 


,5" 


64, 


,8 ±6.9 


1 


.74 ± 





.08 


-1 


.65 


± 


0, 


,05 


J143500+342531 


14h35'"00, 


.98" 


34°25'30, 


,2" 


59, 


,7 ±6.2 


2, 


.18 ± 


0, 


.09 


-1 


.51 


± 


0, 


,05 


J143520+345950 


14h35'"20, 


.51" 


34°59'49. 


,1" 


58, 


,6 ±6.2 


1, 


.55 ± 


0, 


.07 


-1 


.65 


± 


0, 


,05 


J143815+344428 


14h38'"15, 


.28" 


34°44'29. 


.8" 


44, 


,0±4.8 


1, 


.93 ± 


0, 


.09 


-1 


.42 


± 


0, 


,05 


J143331+341012 


14h33'"31 


.84" 


34°10'12. 


.9" 


37, 


,0±4.1 


2, 


.04 ± 


0, 


.09 


-1 


.32 


± 


0, 


,05 


J142631+341557 


14h26'"31, 


.69" 


34°16'00. 


.9" 


22, 


,5 ±3.6 


1, 


.16 ± 


0, 


.07 


-1 


.35 


± 


0, 


,08 


J142954+343516 


14h29'"53, 


.90" 


34°35'18. 


,8" 


16, 


,4 ±2.4 


0, 


.62 ± 





.04 


-1 


.49 


± 


0, 


,07 


J142724+334714 


14h27'"24, 


.77" 


33°47'18, 


,5" 


12, 


,8 ±2.8 


0, 


,43 ± 





.04 


-1, 


.55 


± 


0, 


,11 


J143538+335347 


14h35"^38, 


.76" 


33°53'44. 


,2" 


12, 


,3 ±2.2 


0, 


.56 ± 





.05 


-1 


.41 


± 


0, 


,09 


J143310+333131 


14h33'"10, 


.42" 


33°31'27, 


,7" 


11, 


,3 ±2.2 


0, 


.50 ± 





.04 


-1 


.42 


± 


0, 


,10 


J143230+343449 


14h32'"30, 


.47" 


34°34'49. 


.5" 


10, 


,1±1.8 


0, 


.52 ± 


0, 


.04 


-1 


.35 


± 


0, 


,09 


J142719+352326 


14h27"^19, 


.32" 


35°23'29, 


,2" 


10, 


,0±3.4 


0, 


.42 ± 


0, 


.06 


-1 


.44 


± 


0, 


,17 


J143700+335920 


14h37'"00, 


.74" 


33°59'20. 


.2" 


9.5 ±2.5 


0, 


.42 ± 


0, 


.04 


-1 


.42 


± 


0, 


,13 


J143249+343915 


14'^32'"49, 


.12" 


34°39'14, 


,0" 


6.4 ± 1.8 


0, 


.24 ± 


0, 


.04 


-1 


.49 


± 


0, 


,15 



" From de Vries et al. (2002). 
' Measured at 153 MHz. 



Table 7. 153 MHz catalog selection of 16 USS sources with a spectral index Q153 below —1.3. 



ID 


RA'' 


DEC" 


C) 


1400 (c\ 
•^153 I ) 


^phot 


A 


14h26'"31.75" 


+34°15'57.5" 


-1.48 


-1.35 ±0.08 


4.97 


B 


141^26^47.87" 


+34°58'51.0" 


-0.89 


-0.56 ±0.06 


3.76 


C 


14'^27™41.84" 


+34°23'24.7" 


-0.98 


-0.44 ±0.13 


1.21 


E 


141^32^58.44" 


+34°20'55.4" 


-0.87 


-0.95 ±0.05 


4.65 



"As measured at 1.4 GHz; de Vries et al. (2002). 

* The spectral index between 325 MHz and 1.4 GHz as measured by Croft 
et al. (2008). 

The spectral index between 153 MHz and 1.4 GHz as measured in this 
work. 

Photometric redshift, as measured by Croft et al. (2008). 

Table 8. Spectral indices between 1.4 GHz and 153 MHz for four candidate 
HzRGs by Croft et al. (2008). 



The probability P{i) that candidate i is the true NIR 
counterpart is 



P{^) 



LR,(r, m) 



E,LR^(r,™) + [l-e(mi.J] 



with J running over all possible NIR counterparts. 
is the fraction of radio sources having a NIR counterpart at 
the magnitude limit of the NIR survey. We found 9{m\nn) 
to be 0.62, where miim = 19.2. The values for p{< m) (sur- 
face number density of NIR sources) were estimated from 
the data itself using bins of 0.1 magnitude across the full 
survey area. The a-priori identification fraction 6{< m) was 
also estimated from the data itself, using the technique de- 
scribed by Ciliegi et al. (2003). This involves counting the 
surface number density of NIR sources around the radio 
sources as function of magnitude. This distribution is com- 
pared to a distribution of background objects, the overden- 
sity of NIR sources around radio sources gives an estimate 
of 9{< m). 

We selected 368 radio sources that are present within 
both the 1.4 GHz and 153 MHz catalogs (matched within 
6") and have a simple morphology (fitted with a single gaus- 
sian). We have computed the likelihood ratio of all NIR 



counterparts located within 20" from the radio position. 
We have defined a radio source to have a NIR counterpart 
if LR^ (r, m) > 0.75. The results are shown in Figure 6. 
(3) The identification fraction is ^ 70 percent for > —0.7, 

while for a}|3° < —0.7 the identification fraction drops to 
about 30 percent. This result reproduces the expected (and 
previously observed) correlation between spectral index and 
A"-band identification fraction (e.g., Wieringa 1991). 



5. Conclusions and Future Plans 

We have presented the results from a deep 1.0 mJy 
central noise), high- resolution (26" x 22") radio survey at 
153 MHz, covering the full NOAO Bootes field and beyond. 
This 11.3 square degree survey is amongst the deepest sur- 
veys at this frequency to date. We produced a catalog of 
598 sources detected at 5 times the local noise level, with 
source flux densities ranging from 3.9 Jy down to 5.1 mJy. 
We estimate the catalog to be ^ 70 percent complete at 
the 5.1 mJy flux limit, 92 percent for > 10 mJy sources 
and > 99 percent for > 20 mJy sources, and less than 
1 percent contaminated by false detections. The on-source 
dynamic range (with noise measured near the source) is lim- 
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Figure 6. Identification fraction of radio sources in the near- 
infrared Kg-hand as a function of spectral index between 
1.4 GHz and 153 MHz. 



ited to ^ 600, while ofF-source (noise measured in central 
part of the image, away from bright sources) this rises to 
> 1000. We expect that residual RFI and residual direction- 
dependent calibration errors prevents reaching the thermal 
noise level of 0.2-0.3 mJy beam~^. 

The 153 MHz catalog presented in Section 3.1 allows 
for a detailed study of source populations in a relatively 
unexplored flux range. We have analyzed the source counts 
and spectral index distributions for our survey. From this 
analysis we draw the following conclusions: 

(i) The Euclidean-normalized differential source counts, de- 
termined over a flux range from 10 mJy to 1 Jy, are well 
approximated by a single power-law slope of 0.91. The in- 
consistency with model source counts by Jackson (2005) 
and Wilman et al. (2008) seems due to uncertainty in the 
predicted contribution of FRi galaxies. Our survey is not 
deep enough to detect the flattening at the lowest flux den- 
sities as seen at higher frequencies (Section 4.1). 

(ii) The spectral indices of a subset of 417 sources be- 
tween 153 MHz and 1.4 GHz, limited predominantly by 
the 153 MHz sensitivity, have a median spectral index of 
—0.76. We confirm the flattening trend in the median spec- 
tral index towards lower 153 MHz flux densities, as seen at 
higher frequencies (Section 4.2). 

(iii) The detection fraction of USS sources, having a spec- 
tral index below -1.3, is 3.8 percent, which is comparable 
to the USS fractions of two comparable 153 MHz surveys. 
Six of the 16 USS sources, apparently physically unrelated, 
are found in pairs < 6' apart, which is an unlikely spatial 
distribution to appear at random (Section 4.2). 

(iv) The NIR A'5-band identiflcation fraction of a subset of 
368 sources drops from ^ 0.7 for spectral indices > —0.7 
to ~ 0.3 for spectral indices < —0.7, reproducing a previ- 
ously observed correlation (e.g., Wieringa 1991) which links 
together the known correlations between K-hand magni- 
tude and redshift, and between spectral index and redshift 
(Tielens et al. 1979; Blumenthal & Miley 1979; Willott et al. 
2003; Rocca-Volmerange et al. 2004) (Section 4.3). 

We plan to continue our analysis of the Bootes 153 MHz 
source survey by comparing it with many other available 
catalogs at various spectral bands. A high priority task will 



be to study the properties of the 16 USS sources, to deter- 
mine if these objects are HzRGs, derive estimates of their 
redshifts and search the surrounding area for galaxies at 
similar redshift. This approach has been succesful for the 
identification and study of galaxy cluster formation (e.g., 
Rottgering et al. 1994; Venemans et al. 2002). 

The observations presented here are part of a larger sur- 
vey with six additional, partly overlapping flanking flelds 
observed with GMRT at the same frequency, covering a to- 
tal survey area of ^ 40 square degrees. This same area is 
also covered by extended WSRT observations using 8 bands 
between 115-165 MHz. We will combine these observations 
with the observations presented here to produce a com- 
bined high- and low-resolution catalog at ^ 153 MHz to 
further facilitate the study of the low-frequency sky, and 
in particular to facilitate the further search for USS radio 
sources. 
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Appendix A: Selection of 153 MHz Radio Source Images 



FLUX = 543.9rnIy F(MS = 2.aBmIy/beai 



FLUX = 335.9 inJy RMS= 1.24 mJy/bs; 



iJy (iMS = 1.18rnly/beam 



FLUX = 149.4 miy RMS = 1.08 mJy/bearr 




FLUX = 32.9 mJy RMS - 1 .OB nUy/bi 
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Figure A. 1. Grey-scale images of a selection of resolved 153 MHz sources with a complex morphology (i.e., not a point source or 
double point source morphology). Horizontal and vertical axes are RA and DEC, respectively. For each source, the total flux and 
local RMS noise are specified above the 5' x 5' image. Contours are drawn at [—3, 3, 5, 10, 20, 50, 100] times the local noise. The 
grey-scale ranges from —3 to 30 times the local noise. 



